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Abstract 

The relation between algebraic and traditional calculations of molecular vibrations is investigated. 
An explicit connection between interactions in configuration space and the corresponding algebraic 
interactions is established. 
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Ab initio calculations for rovibrational spectra of molecular systems attempt exact solutions of the 
Schrodinger equation. In practice, the molecular Hamiltonian is usually parametrized as a function of 
internal coordinates |Q and the potential is modeled in terms of force- field constants, which are determined 
through calculations involving several configurations associated to the molecular electronic states g] . For 
small molecules this procedure is still feasible, but this is in general not the case for polyatomic molecules, 
due to the large size of the configuration space. It is thus important to develop alternative methods to 
describe these systems. Algebraic (or vibron) models attempt to provide such alternative techniques || [|] . 
In its original formulation [|5], ^| rotations and vibrations were treated simultaneously in terms of coupled 
£7(4) algebras: Q = t/i(4) ® ^(4) 8> ■ • • • For polyatomic molecules it was found to be more convenient to 
first separate the rotations and vibrations and subsequently to treat the vibrations in terms of coupled 
U(2) algebras: Q = U\(2) ® £^2(2) <£> . . . . In the latter version of the vibron model the calculation of 
matrix elements is greatly simplified. An additional advantage is that it is well-suited to incorporate 
the underlying discrete symmetries H, |[ [|. The vibron Hamiltonian, however, is expressed in terms of 
abstract algebraic operators, whose connection with more traditional methods has been limited to studies 
of the corresponding energy surface |M . 

The aim of this letter is to investigate the relation between algebraic and configuration space calcula- 
tions and to establish an explicit connection between interaction terms in coordinate space parametriza- 
tions and in the more abstract algebraic space. We shall illustrate this connection by studying the Be4 
cluster, for which the force field parameters were determined in [ p"l| by ab initio methods. Our compari- 
son is based on an analysis of tetrahedral molecules in terms of symmetry-adapted internal coordinates 

©■ 

In the algebraic approach each relevant interatomic interaction is associated with a U(2) algebra. 
In the present example of the Be4 cluster, which has a tetrahedral shape, there are six U(2) algebras 
involved: Q = Ut(2) ® U 2 {2) ® U 6 {2). Each t/,(2) algebra (i = 1, ... ,6) is generated by the set 

{Gi} = {Ni, J+,i, J-,i, Jo,i}, satisfying the commutation relations 

[Jo,iiJ±,i] — ±J±,i 5 [J+,i,J-,i] = 2Jo,i , [iVj, J^i] =0, (1) 

with fi = ±, 0. Here Ni is the (boson) number operator and the operators J M .; satisfy the 'angular 
momentum' commutation relations of SUi{2). Since Jf — Ni(Ni + 2)/4 we can make the identification 
ji = Ni/2. The eigenvalues of Jo,i are restricted to m, > and can have the values mi — Ni/2, (Ni — 
2)/2, . . . , 1/2 or for Ni odd or even, respectively |0]. The local basis states for each oscillator are usually 
written as \Ni,Vi), where Vi = (Ni — 2rrii)/2 = 0, 1, . . . [Ni/2] denotes the number of oscillator quanta in 
the i-th oscillator. Because of the tetrahedral symmetry of the Be4 cluster Ni — N for the six oscillators. 
The operators in the model are expressed in terms of the generators of these algebras, and the symmetry 
requirements of the Hamiltonian under the tetrahedral group Td can be readily imposed [D| . In the 
usual algebraic formulation different chains of subgroups of Q are considered and the Hamiltonian is built 
by means of appropriate combinations of invariant operators associated to these chains 0, ^ ^) ■ Here 
we shall follow a different approach which leads to a richer structure for the Hamiltonian and to a direct 
connection to configuration space interactions. To achieve this goal we note that the local operators {Gi} 
acting on bond i can be projected to any of the 7^ fundamental irreps T = A\, E and Fi. Using the J Mi i 
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generators of Eq. (|l]) we obtain the Td tensors 

6 

^£,7 = X! "7^ ^ ' ( 2 ) 
i=l 

where /x = ±, and 7 denotes the component of T. The explicit expressions are given by 
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1 

2^3 



t=i 



- -j=[Jn,3- Jv.,5) ■ (3) 



V2 
1 

72 

The algebraic Hamiltonian can now be constructed by repeated couplings of these tensors to a total 
symmetry Ax, since it must commute with all operations in Ti- This is accomplished by means of the Ti 
Clebsch-Gordan coefficients ||, 

In order to establish a connection with configuration space calculations we use the analysis of 
In this work the vibrational Hamiltonian for the Be4 cluster is expressed in terms of symmetry- adapted 
internal coordinates, q^, and momenta, p^. The transformation to the tensor operators of Eqs. (^||) 
proceeds in two steps. First we introduce creation and annihilation operators 

Here the normal boson operators are related to the local boson operators by b^ — J2i al yi^i- The 
coefficients o£. i can be read from Eqs. (||||). Next the local boson operators are associated with the 
generators of Eq. (|l]) by means of 

h -» J+,i/VNi , b\ -> J-,i/y/Ni ■ (5) 

This transformation is such that the commutator 

j^^J-A = ±*So,i = = l-f , (6) 

reduces for iVj — ► 00 to the standard boson commutator &]] = 1. Eq. (^|) provides a procedure 
to construct an anharmonic representation of harmonic operators. The anharmonic contribution arises 
from the —2%i/Ni term in Eq. Note that each local harmonic oscillator (b\bi + bib\)/2 leads to 

J+,j + J +: iJ-s)/2Ni — Vi + 1/2 — vf/Ni and thus to a Morse-like spectrum through the association 
of Eq. (||). In this way the algebraic model substitutes harmonic oscillators by anharmonic ones. 

This procedure can be applied to the various contributions to the vibrational Hamiltonian for the Be4 
cluster |0|. The vibrational basis states for this system are usually labeled by (fx, z/™, v l 3 ) fl4j] . Here v\, 
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i>2 and z/3 denote the number of phonons in the A\, E and Fi modes, respectively, and m = 1*2, v-i — 2, . . . , 1 
(or 0) for vi odd (or even) and and I — z/3, 1/3 — 2, . . . , 1 (or 0) for 1/3 odd (or even) are the vibrational 
angular momenta associated with the E and F% modes. For the zeroth order vibrational Hamiltonian we 
find 

H = uji H Al +u 2 H E + ^3 7-Lf 2 i (?) 

with 

"r - ^E( f -,^ + n,^) ■ (8) 

7 

The anharmonic vibrational terms are expressed in terms of products of TLy > 

Hi = Xn (Hai^ + X22 (He^ + X33 [Hf^ 

+X12 {Hax'Hej + -X"i3 ( 'H-A^'H-Fz^j + X23 (HeHf 2 j ■ (9) 

The further splitting of vibrational levels (^1,^2,^3) into its possible sublevels is achieved by means of 
the interactions |0] 

7 

+*33 [ 6^[f F2 x T F *] F [Tf x ] F - 4^[T Fa x f F2 ]f [Tf 2 x f f}^ j ^ 

V 7 7 / 

+i 23 (8^[f F x f F2 ] Fl [if x f F2 ] F ^ - 8^[f F x T F ^ [ff x f?]A JL 
\ 7 7 / 



(10) 



The operators and l Fl represent the vibrational angular momentum operators associated with the E 
and F 2 modes, respectively, 



l M = -iV2-[f F x ff] A2 , 



f Fl = +i\/2l[f F2 xf F2 ] Fl . (11) 
The square brackets in Eqs. ([Io|,[ll]) denote the tensor coupling under the point group Td 

[f T ^f T % = J2 ^(ri,r 2 ,r ;7l , 72 , 7 )f 7 r il f 7 r 2 2 , (12) 

71.72 

where the expansion coefficients are the Clebsch-Gordan coefficients for Ti || [L3). The interactions of 
Eq. ( p^ ) were absent in previous versions of the model § §, §. 

The algebraic Hamiltonian of Eqs (|7|-|l2]) is the algebraic equivalent of the vibrational Hamiltonian 
of p2[ . The harmonic frequencies u>i and anharmonic constants Xy, 322, 533, £33 and t 2 3 have the 
same meaning as in jl2) and the ab initio calculations of jllj can be used to generate the spectrum. 
The various contributions to the algebraic Hamiltonian arise naturally from the successive couplings of 
the fundamental tensors of Eq. (||). The scale transformation of Eq. (||) makes it possible to establish 
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the connection between ab initio and algebraic parameteres and to explicitly construct the algebraic 
interactions that correspond to interactions in configuration space. In the opposite sense, Eqs. (j^^) 
provide a procedure to obtain a geometric interpretation of algebraic interactions in terms of those in 
configuration space. In the harmonic limit, which is defined as iVj — > oo, Eq. (Q) reduces to the standard 
boson commutator [E>^ ? &| ] = 1. This limit corresponds to a contraction of S77j(2) to the Weyl algebra, 
generated by the set In the harmonic limit the Hamiltonian Hq + Hi + H 2 of Eqs. |l2|) 

reduces exactly to the vibrational Hamiltonian of |Q . Note that because of the replacement of Eq. (||) 
the Hamiltonian Hq + Hi + H2 only depends on the T^ 7 tensors of Eq. (|^) with fi = ±. In addition, 
the algebraic model provides terms involving the Tj* tensors with /1 = 0. As can be seen from Eq. (||) 
these terms are completely anharmonic in origin and have no direct counterpart in models based on the 
standard harmonic bosons. These operators arise from the substitution of harmonic oscillators by Morse 
oscillators and play an important role when dealing with anharmonic molecules, particularly at higher 
phonon numbers [fj"5| . 

Apart from providing a direct connection to configuration space calculations, this formalism can also 
be used as an effective model of molecular vibrations, particularly when no ab initio calculations are 
available. As an example, we show in Table | the results of a fit to the ab initio calculations for Be4 up 
to four phonons. The ab initio results were generated with the parameters from |ll| . The Hamiltonian 
used in the fit contains 9 interaction terms compared to the 13 of Q (see Eqs. (|7|,|9|,|l0|)). The parameters 
are extracted in a fit that includes all vibrational energies up to four phonons (V = v\ + v 2 + v 3 < 4): 
u) X = 636, lo 2 = 453, to 3 = 53 2, X 33 = 44.2 76, X 12 = 4.546, X 13 = -2.539, 533 = -15.031, t 33 = -1.679 
and £23 = —1.175. All values are given in cm' 1 . The total number of bosons used in the fit is N = 44. 
The r.m.s. deviation between the algebraic and the ab initio calculations is 2.6 cm -1 . 

As a test of the predictive power of the algebraic approach we have performed another calculation 
in which the same 9 parameters were determined in a fit that only included the vibrational energies up 
to three phonons (V < 3). In this case the r.m.s. deviation is 1.6 cm -1 . If we now use these values of 
the parameters to calculate the four phonon states, the r.m.s. increases to only 3.1 cm -1 , compared to 
2.6 cm -1 in the previous calculation. We remark that by restricting the model interactions to Casimir 
invariants and their powers J|, ||, f?J the Be4 spectrum cannot be reproduced. The terms in H 2 of Eq. ( |io| ) 
play a crucial role. 

Repeating the same fit in the harmonic limit [N — ► 00) the r.m.s. deviation increases from 2.6 to 5.4 
cm -1 . This shows that the anharmonic contributions introduced by taking a finite value of N (see Eq. (^|)) 
provide an important improvement of the fit. The real test of this aspect is a fit to experimental data 
rather than to other calculations. Work on the application of this algebraic model to the experimental 
vibrational spectra of polyatomic molecules is in progress Q . 

In summary, in this letter we have established a connection between algebraic and configuration-space 
interactions. For the example of the Be4 cluster (with tetrahcdral symmetry) we have, starting from 
configuration space interactions, constructed explicitly the corresponding algebraic interactions (which 
have a richer structure than in previous versions of the model). In the harmonic limit the configuration 
space results are reproduced exactly. 

In addition, it was shown that the algebraic model can also be used as an effective model of molec- 
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ular vibrations with good precision. In the algebraic approach the eigenvalues and corresponding wave 
functions are obtained by matrix diagonalization. Hence the required computing time is small. These 
properties open the possibility to use the algebraic model as a numerically efficient, empirical tool to 
study molecular vibrations, especially when no ab initio calculations are available (or feasible) fll5|] . 

We thank F. Iachello and P. van Isacker for their continuous interest and useful comments. This work 
was supported in part by the European Community under contract nr. CIl*-CT94-0072, DGAPA-UNAM 
under project IN105194, CONACyT-Mexico under project 400340-5-3401E and Spanish DGCYT under 
project PB92-0663. 
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Table I: Vibrational excitations of Be4 using the algebraic Hamiltonian with parameters given in the text. The 
ab initio (JV — > oo) spectrum is generated with the parameters from JL1|. The energies are given in cm -1 . 
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Tabic I: Continued. 
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